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1. Introduction 

The dynamics of interacting agents is currently studied intensively, applying the ideas 
and techniques of equilibrium and non-equilibrium statistical mechanics. One of the 
models which has attracted particular attention is the so-called minority game (MG), 
introduced as a minimalist econophysics model for a financial market pQ. The players 
in the MG are traders who, at each round of the game, have to make one of two 
possible choices (e.g. buy or sell) in response to publicly available information. Each 
aims to make profit by making the opposite choice to the majority of agents. The 
interaction between agents is indirect: they cannot observe individual actions of others, 
but only the subsequent cumulative effect of all actions on the market. To determine 
their own trading actions, each agent holds a pool of strategies, assigned randomly 
before the start of the game and then kept fixed. These effectively act as look-up tables, 
mapping the observation of publicly available information onto a proposed trading 
action. In the versions of the MG studied so far (e.g. PQ 121 EJ E]) agents cannot 
combine strategies, but select the one which they regard as their best. We refer to those 
types of MGs as conventional. The identification of the best strategy is based on points 
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the agents allocate to each of their strategies in order to measure their performance. 
After each round of the game each agent evaluates the quality of each of his or her 
strategies, increasing the points of strategies which would have yielded a correct minority 
prediction. For a general overview of MG-type models we refer to [3] . 

The update rules of the MG look simple, but describe surprisingly complex 
cooperative processes. This is most visible in the non-trivial behaviour of the 
fluctuations of the total bid, the so-called market volatility [TJ 13 12 El- The main 
control parameter in MGs is the ratio a = p/N of the number p of possible values for 
the public information over the number N of players. One observes a critical value a c 
which marks a dynamical^ phase transition, separating a non-ergodic phase (a < a c ) 
from an ergodic one (a > a c ). In the non-ergodic phase, the volatility is very sensitive 
to initial conditions [211112], and the integrated response is infinite j3j. Moreover, in 
the stationary state the system exhibits persistent oscillations in the non-ergodic phase, 
whereas oscillations decay on finite time scales for a > a c jHCHl- 

Analytical progress is possible using equilibrium and dynamical approaches and has 
resulted in analytical expressions for a c , which are now regarded as exact [TT) 1^ IT^ IT3] . 
The generating functional analysis a la De Dominicis ^3] has proven particularly 
valuable; it enabled a full understanding of the dynamics of the MG in the ergodic 
phase. In this formalism, the strategy selection dynamics of the agents is mapped 
onto a non-Markovian effective single agent process. In the case of conventional MGs, 
the microscopic laws and the resulting single-trader process are non-linear and resist 
analytical solution. Instead, one derives a coupled set of implicit equations for stationary 
states, from which one tries to extract the values of the persistent order parameters. In 
contrast to equilibrium systems, there are no fluctuation-dissipation relations which 
could be used to simplify those equations. In the ergodic phase the analysis can 
be simplified taking into account the existence of so-called 'frozen agents' (runaway 
solutions of the microscopic laws). A proper understanding of the dynamics in the non- 
ergodic regime, however, is still lacking. Moreover, for conventional MGs the market 
volatility (the MG's main observable) cannot be expressed in terms of persistent order 
parameters. Instead, detailed knowledge of both the long-time and the short-time 
behaviour of the macroscopic order parameters is required. Hence, even in the ergodic 
phase, results for the volatility are so far restricted to approximations, whereas in the 
non-ergodic phase only approximate asymptotic results in the limit a — > are available 
[3]. For a recent review of dynamical MG analyses see e.g. ^Hj- Approximations for the 
volatility in the ergodic state are also accessible within the framework of replica theory 

In this paper we present a version of the MG which is analytically solvable, but 
nevertheless displays some of the interesting features found in the conventional MG. To 
this end we study the dynamics of a spherical version of the MG using the generating 
functional approach. Like in spherical p-spin glasses with polynomial equations of 

| MGs do not obey detailed balance, so one can only speak about non-equilibrium phase transitions. 
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motion for the continuous microscopic degrees of freedom, explicit closed equations for 
the two-time correlation and response functions can be formulated |16| . 

A second control parameter, the radius r of the sphere to which the dynamics is 
confined, becomes relevant in the present model. Apart from the spherical constraint 
we choose the update rules to be linear in this paper, so that we can solve the resulting 
dynamical equations exactly, reminiscient of the p = 2 case known for spherical spin- 
glasses |17j . In particular we are able to compute the volatility in all regions of the 
phase diagram without making any approximations at any stage. In terms of the decision 
making of the individual agents the linear spherical model corresponds to allowing them 
to play linear combinations of strategies (rather than to pick the best one). 

Despite its simple microscopic rules, the spherical MG as presented in this paper 
shows interesting behaviour and exhibits novel features as well as properties analogous 
to the ones of conventional MGs. In particular we find three distinct phases in the 
(a,r)-plane. Our analytical findings are verified convincingly by numerical simulations. 

2. Model Definitions 

Before defining the spherical version of the MG and giving an interpretation of its update 
rules, we will recall the dynamical rules of a conventional MG as studied for example 
in [IJ ^]. We label the N agents in the MG with Roman indices. At each round t 
of the game each agent i takes a trading decision 6j(t) G IR (a 'bid') in response to 
the observation of public information 1^ which is chosen randomly and independently 
from a set with p = aN possible values§, so //(£) G {1, . . . ,aN}. The rescaled total 
market bid at round t is defined as A(t) = N~ 1 / 2 ^ h(t). Each agent i has S > 2 
fixed trading strategies (look-up tables) Rj a = (R} a , ■ ■ ■ , Rf a r ) at his or her disposal, 
with a — 1, . . . , S. If agent i decides to use strategy a in round t of the game, his or 
her bid at this stage will be bi(t) = Rf^ ■ All strategies Rj a are chosen randomly and 
independently before the start of the game; they represent the quenched disorder of this 
problem. The behaviour of the MG was found not to depend much on the value of S 
fHlEU], nor on whether bids are discrete or continuous [2|. For convenience, we choose 
5 = 2 and H ia G { — 1, l} aAr in this paper. In order to decide which strategy to use, the 
agents assign points Pi a (t) to each of their strategies, on the basis of what would have 
happened if they had played that particular strategy: 

Pia(t)^p m (t)-R^A(t). (1) 

Strategies which would have produced a minority decision are thus rewarded. In the 
conventional MG, at each round t each player i uses the strategy in his or her arsenal 
with the highest score, i.e. b^t) = R^~%y where a^t) = arg max a p ia (t). For 5 = 2 the 
rules ((H) can then be simplified upon introducing the differences qiit) = \[pnit) — Pi2(t)}. 
Thus, if qi(t) > 0, agent i plays strategy R^, whereas for qiit) < he or she plays R«2- 

§ This is the so-called MG with random external information; in the early MG definition JS| the 
external information was not random but coded for the actual history of the global market. 
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Hence, in the conventional MG bi{t) = cuf + sgn[(fe(t)]£f , where u)i = §[Ra + R*2] 
and £ 4 = |[Rji — Rj2]- The evolution of the {g^} is given by 



m(*) 



(2) 



with O = N~ 1 / 2 Y2j Uj. Equation (J2J) defines the standard (or so-called 'on-line') MG. 
Alternatively, corresponding to updating the {q{\ only every O(N) time-steps, one might 
define the dynamics in terms of an average over all possible values of the external 
information in (J2J), resulting in the so-called (conventional) 'batch' MG 

qi(t + 1) = qiit) Jij sgn[gj(t)]. (3) 

3 

Here Jy = 2N~ 1 £ i ■ £ ■ and /ij = 2iV _1 / 2 £j • f2. See [T3] for stochastic extensions and 
[HI] for consideration of the effects of anti-correlation of strategies on the comparison of 
on-line and batch models. 

The batch model (jHJ) is particularly suitable for being replaced by a spherical 
version. We first linearize Q, and subsequently normalize)! the vector (qi, . . . ,q^) 
to a fixed length r > at each iteraction step t, resulting in the spherical batch MG: 

[1 + X(t + l)] qi (t + 1) = qi (t) J ij Qj(t), ( 4 ) 

^E^) = r2 foralH - ( 5 ) 

i 

The values of the constraining forces X(t) in (J3J) follow from (0); we exclude artificial 
sign changes by insisting on 1 + \(t) > for all t. We note that our model ()4I5|) has no 
analogue of the tabula rasa MG initialization, qi(0) = for all i, often employed in the 
conventional MG. We also note that the {q{\ of the conventional MG (0) do not satisfy 
a spherical constraint, unlike the spins of a conventional Ising spin system. There is 
thus no reason for restricting oneself to r = 1. In fact, r is a new control parameter 
and the system exhibits phase behaviour in the (a, r)-plane with interesting differences 
from the conventional game. 

The linearity of (JIJ implies that agents now play linear combinations of their 
strategies. Upon presentation of public information at time t the bid of player % 
in a corresponding on-line game is 

hit) = i[i + qMK + i[i - ft(*)M- (6) 

The main object of natural interest in MGs is the volatility, which describes the standard 
deviation of the total (re-scaled) market bid 

^[q(*)] = ^E^ + ®W]- ( ? ) 

* i 

|| The spherical normalisation is necessary to suppress possible runaway solutions corresponding to 
cigcnmodcs of the linear update rule with eigenvalues of a modulus larger than one. 
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In the on-line models the relevant averages are over the stochasticity of the 'information'. 
In deterministic batch problems, such as discussed here, these averages are replaced by 

ones over //: (At) = P~ x EJ=i ^IqW] and (A t A t >) = p~ x E£=i a *W)\ a *W)\- The 

jl] and define a more general 



volatility is defined as of = (A 2 ) 
object, the volatility matrix S tt / = 



- (A t ) 2 . Here we follow 
(A t A t ,)-(A t )(A t ,): 



p P P 

- A"[q(f)]A"[q(f )] - \ E A"[q(t)] - AW)) 



(8) 



Note that of 



E#. Random trading, with q(£) taken randomly and independently 
from the sphere q 2 (t) = Nr 2 at each time t, would result in (A t ) = 0(N~^) and 
Et# = | + \ r2 $tt' + 0(iV~^). The volatility measures the efficiency of the market, with 



07 = corresponding to a perfect match between supply and demand at time t. 



3. Macroscopic Dynamics 

The similarity between the spherical batch MG ()4I5|) and the conventional batch MG 
(jSJ) allows us to obtain the effective single trader equations for ()4I5|) simply by making 
the substitutions q(t + 1) — > [1 + A(£ + l)]q(t + 1) and sgn[g(t)] — > g(t) in the results of 
0] (found within the generating functional analysis framework, in the limit iV — > 00): 

[1 + \{t + l)]q(t + 1) = q(t) + 9{t) - a ^(1 + G) u }q(t') + ^ T]{t). (9) 

t' 

Here is an external perturbation field introduced to generate response functions and 
f](t) is a zero-average Gaussian noise, characterized by the following covariance matrix 
(with D tt , = 1 + C tt , and l tt ' = 8 W ): 

Z tt , = (v(t)v(t% = [(l + G)- l D(l + G T )- l ] al . (10) 

The matrices C and G and the constraining forces X(t) are the dynamical order 
parameters of the problem, to be determined self-consistently by solving 

C tt > = (q(t)q(f))., Gw = -^-(q{t))^ C tt = r 2 . (11) 

One always has (q(t))* = 0. The physical meaning of C and G is given by 

Jv— >oo iv ' 

i 

i 

where 7TT denotes an average over the disorder, i.e. over the space of all strategies in 
the context of the MG. The brackets (...)* in |T0|) and fjl 1 jl refer to averaging over the 
realisations of the process (JOJ), i.e. over the noise {rj{t)}. As usual, the single-agent 
process © is non-Markovian. 

It is possible to convert the system ()9Jl()Jllj) into a pair of explicit iterative 
equations for C and G. An explicit equation for C results upon multiplying Q 
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by q{t') and subsequently averaging over the noise. We make use of the identity 
(r](t)q(t'))z = V^EsPA) (derived via an integration by parts in the generating 
functional; see [3] for an analogous identity in the conventional MG). To deal with C 
(which obeys causality: G tt > = for t < t') one takes a field derivative of (JHJ), followed 
by averaging. The result reads: 

[1 + \{t + 1)] Ct+i )t , = Cut +a[(l + G)- 1 J D(1 + G T y 1 G T ] tt ' 

-a[(l + G)- 1 C] tt ,, (14) 
[1 + X(t + 1)] G w > = G tt , - a[(I + G)~ l G\ tt > + 6 W . (15) 

These coupled equations have to be solved subject to the constraint C tt = r 2 for all 
t > 0. Furthermore, as in [4j one finds for N oo that the rescaled disorder-averaged 
average bid (At) is zero at any time, and that the disorder-averaged volatility matrix 
(JHJ) is proportional to the covariance matrix of the single-trader noise: 

lim 5* = ±53* = \ [(l + Gy'Dil + CFyX,. (16) 

The dynamic order parameters in the spherical model are prescribed in full by 
([HUl]) , wr th the constraints C u = r 2 and 1 + A(t) > 0. As in jl], they can be calculated 
iteratively, starting from (t,t') = (0,0), and upon using causality (i.e. [G n ]^/ = for 
n > t — t') . Given the prescribed values Coo — r2 an d Goo — one finds, for instance: 

A(l) = - 1 + (1 + a(r~ 2 - 1) + a 2 ) 1 / 2 , (17) 

Gjo = [l+a(r- 2 -l) + a 2 ]"^, (18) 
r 2 (l - a) 

C 10 = — 1 ' (19) 
v/l + a(r- 2 - l) + a 2 

(note: G\\ = and Cu = r 2 ). From these follow the volatility matrix elements: 

Soo = 1 + r 2 , (20) 
ar 2 + 1 

Sio = 1 - , (21) 

^l + a(r~ 2 -l) + a 2 

? 2 l-r 2 (l-2a) , , 

VJ _ -y _|_ r ^ I V I (22) 

11 " ^l + a{r~ 2 - l) + a 2 1 + a(r- 2 - 1) + a 2 ' 

This iteration can be carried out for an arbitrary number of time steps; in practice, 
however, the terms become prohibitively more complicated for increasing times. 
Alternatively, one may iterate (|14|15|l numerically. Unlike the procedure proposed in 
[20] to generate realisations of the single trader process, the iteration of (|14|15jl does 
not require averaging over the single agent noise, and thus provides very precise data 
(albeit that numerically inverting the t x t matrices in ()14|15j) becomes more and more 
costly as the number of time steps increases). 

One observes that, due to the explicit form of ()14|15|) (and in sharp contrast 
to the similar calculation in jlj), the temporal evolution of the macroscopic order 
parameters is completely independent of initial conditions: as long as the constraint 
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liniTv^oo A -1 Yli=i <?i(0) 2 = Coo — f 2 is me t the distribution P(gj(0)) from which the 
initial point differences are drawn is irrelevant for the values of both C and G at finite 
times as well as for the macroscopic stationary state. We have verified this in numerical 
simulations, initializing the dynamics with different distributions for the gi(0), but all 
with second moment A -1 ^ i=1 qf(Q) = r 2 . This property of the spherical MG is quite 
distinct from the conventional MG, where the explicit analysis of the first few time 
steps as presented in jl] reveals that the values of the correlation and response functions 
at finite times depend on the higher moments of P(gj(0)) as well, and not only on its 
variance. As far as the stationary state of the conventional MG is concerned the interest 
so far has mainly focussed on starts of the form \qi(0) \ = go for all i, with go — for tabula 
rasa starts and go > for biased starts. Crucial differences between the two cases have 
been found in the non-ergodic regime [HI El EH- We have extended this analysis and have 
performed simulations of the conventional batch game with different initial distributions 
for the <7i(0) all with the same second moment q 2 . For fixed second moment we find 
that qualitatively the stationary volatility does not depend on the higher moments of 
P(gi(0)), but that differences in the quantitative values are found for different shapes of 
the initial distribution of point differences. 

Let us finally inspect the solution of (|14I15|) for small a and finite times. By 
induction one finds 

X(t) = 0(a), G tt > = &w + 0(a), C tt> = r 2 + 0(a), 

where Q t f = 1 if t > t' and Q tt i = otherwise. For the volatility we find E tt = O (a 2 ) 
for all finite t > 2. We conclude that for small a the system is in a completely frozen 
state, with a divergent integrated response and vanishing volatility. 



4. Analysis of Stationary States 

4-1. Implications of Time- Translation Invariance 

We now focus on time-translation invariant solutions of the dynamical equations (|14|15j) . 
i.e. we consider the system long after any initial equilibration and study solutions of 
the form 

C w = C{t-t), G tt >=G(t-t'), \{t)=\. (23) 

Then all matrices in (jl4H5j) become Toeplitz matrices, and hence they commute. Given 
(j2HJ) it is natural to express the dynamics in terms of the Fourier transforms of the 
correlation and response functions. We use the following notation: 

C(t) = f ^e^C(co), C{u) = e-*"C{T) (24) 
and similarly for G. The equations ()14|15j) subsequently translate into 

k{ U )C{ U ) = a »M?W - ^ , (25) 
V ; V ; |l + G(u;)| 2 1 + G(u) 
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Figure 1. Phase diagram of the spherical MG, displaying three phases: O, 
oscillating correlation function and finite integrated response; F, frozen phase with 
finite integrated response; and AF, anomalous frozen phase with diverging integrated 
response. Throughout the phase AF the volatility is zero. The O— >F transition is 
defined by (|39[) . and the F— >AF transition is defined by 146(1 : both are continuous. 
The discontinuous transition from O— >AF occurs at a = 3 — 2^/2 w 0.172. The triple 
point corresponds to a = 3 — 2\/2 and r — r* « 0.455. 



AMGM = 1 - ^H- (26) 
1 + G[oj) 

where A(lj) = (1 + X)e lw — 1, and G(u)* denotes the complex conjugate of G{u). 
Since D{ui) = C{oj) + 2w8(uj), and upon defining the integrated static response 
X = J2 T ^( r ) = ^(0)) we ma y rewrite (|23j) as 

A(w)|l + G(u)\ 2 + a\c{u) =2ira X 5(u). (27) 

Note that considering the case u = in (}26|) allows us to express A in terms of x- 



A = 1+ * (1 ~ a) . (28) 

x(i + x) 

Finally, in the stationary state the volatility matrix 5 is also of the Toeplitz form 
E u i —E(t — t') and thus the volatility a 2 = S(0) can be expressed as 

2 1 f n duj C(u) 1 



a 2 J /_ 7r 2 7 r|l + GMP + 2(l + X )2- (29) 

Since initial conditions play no role in the macroscopic dynamics, we must conclude that 
as soon as multiple stationary solutions exist only one of these will ever be realized. 



Dynamics of a spherical minority game 



9 



Before we give a detailed account of the further analysis of the dynamical equations, 
we briefly summarise our results. We find that, depending on the control parameters a 
and r, the system displays three distinct phases, as illustrated in figure HJ 

(i) a phase with finite integrated response and oscillatory behaviour of the correlation 
function (O), 

(ii) a frozen phase with finite integrated response (F), 

(iii) a frozen phase exhibiting anomalous integrated response which grows linearly with 
time (AF). 

We will now proceed to obtain exact solutions of the dynamical equations in each of the 
three phases. 

4-2. Stationary States without Anomalous Response 

We first inspect stationary states with finite Xi i- e - those for which perturbations will 
decay sufficiently fast. It now follows from (|2?j) that, for any uj 7^ 0, C(oj) can be non- 
zero only if A(a;)|l + G(u)\ 2 + a = 0. This requires A(a;) to be real, which is possible 
only for u = 0, n. We conclude that C(u) = 2ttc 5(uj) + 2ixci5{uj — 7r), or equivalently 

C(r)=co + ci(-l) T . (30) 

Here Cq and C\ (which will depend on the parameters a and r, as will x) are coupled via 
the spherical constraint: Cq + C\ = r 2 . Insertion of (J3*U|) into f)27j) leads to the following 
two coupled equations 

c [a + A(l + x) 2 } = ocx, (r 2 ~ c ) [a - (A + 2)(1 + X ') 2 } = (31) 

with x' = — = C( 7r ) measuring the response to persistent oscillating 
perturbations. These equations are to be solved in combination with We conclude 
that there are two types of stationary states with finite x : a frozen state, where Cq = r 2 
(so Ci = 0), and an oscillating state, where cq < r. We will work out their properties 
separately below. For solutions of the form (JHUj) we can also work out (J29~j) further: 

2 1 +C (r 2 - Cq) 

= WM , no + o/1 , -.AS ' ( 32 ) 



2(1 + X ) 2 2(1 + X ') 2 ' 

Oscillating stationary states without anomalous response. Here Cq < r 2 , and the 
remaining four (coupled but closed) equations to be solved to find the stationary state 
include one expression for x' which one obtains by choosing u = n in (|26j) : 

<^X r\ , oui , ,.A2 



co= - , "r, ,, 2 , a = (A + 2)(l + X ') 2 , (33) 



, l + y(l-a) , l + x'(l-a) , s 

The set f!33l34j) allows for two types of solutions. The first, where A = a — 1 + 2^/a, 
obeys the requirement 1 + A > for all a. One must in fact demand A > in order to 
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have a finite x ( as required), and we reject x < solutions on physical grounds. This 
leaves: 

A = a-l + 2^, (35) 

l-a-y/a + V2a 3 / 2 + a 2 

X = + a ' (36) 

1 + \/at 



a c ,2(r) 



(39) 



We note that A, x, x' an d Co are independent of r; only cj depends on r via ci = r 2 — Co. 
The second type of solution, where A = a — 1 — 2y/a, meets our requirement A + 1 > 
only for a > 4. It turns out that such solutions are never realized, so we will not give 
their equations in full. We have now determined all order parameters and the volatility 
in explicit form: cq follows from insertion of ()35I36|) into the first equation of (}3*3*j) . 
whereas the volatility follows upon inserting Co and (|3fi|37jl into ((HZj) . For a — > oo one 
finds lim^oo x = 0, lim^oo A/a = 1, so that limQ^ooCo = 0; the amplitude of the 
oscillations in the correlations increases with increasing a. 

The present solution breaks down when either x oo or C\ — > 0. The corresponding 
mathematical conditions are found to be a = a Cj i (with x < oo for a > a c ,i) and 
ot = a C) 2(r) (with c\ > for a > a c ^(r)), respectively, where 

a cA =3-2V2 w 0.172, (38) 

2 + 1/r 2 
2y/l + l/r\ 

We note that a c ,i = a C) 2(r) at r = r* = \Ja c ,\/ (1 — tt c ,i) ~ 0.455, and that a Cj2 (r) > a c>1 
for r < r* and a C) 2(r) < a Cj i for r > r*. Thus one expects that, as a is lowered for any 
r > r*, the amplitude of the oscillations of the correlations remains positive until the 
critical value a Cj i is reached and a transition to a state with anomalous response occurs. 
At this point one has 

lim a 2 = 1 — = ( r 2 - Jj^L ] . ( 40 ) 

aiac.1 3-2^^ 2(^-1)/ 

For r < r* the oscillatory behaviour of the correlations breaks down as a is lowered 
before anomalous response sets in, and the system enters a frozen state with finite 
integrated response. 

Frozen stationary states without anomalous response. Here Cq = r 2 , and there is no 
need to calculate x'i the coupled equations to be solved are simply 

Ax 2 + (2A-a/r 2 )x + A + a = 0, (41) 
Ax 2 + (a - 1 + X) X - 1 = 0, (42) 
from which we can determine A and x to be: 
(2 + l/^)V5 
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x = [vW* + V^ 2 - i]" 1 - (44) 



Upon inserting (jUjl and the relation c = r 2 , one finds that our exact expression 
for the volatility simplifies to 



1 r i 2 



Vr 2 + 1 - r/v^ . (45) 

The present frozen state will cease to be a consistent solution at the point where \ 
diverges. Equation (}4*4"|) states that this happens at a = a c ^(r), where 

a c ,s(r) =r 2 /(r 2 + l). (46) 

The line a = a Cj 3(r) marks the transition between a frozen state with finite integrated 
response (for a > a c ,3 (?")), and an anomalous frozen state (for a < 0^3 (r)). According 
to (|45)l it also coincides with the line where the volatility vanishes. 



4-3. Frozen States with Anomalous Response 

Numerical simulations for small a and large r reveal a parameter regime with a 
stationary state in which the volatility vanishes, a 2 = 0, and where all agents are 
frozen in such a way that 

C(r)=r 2 , A(r) = 0, X = oo. (47) 

These three identities clearly hold at the transition line (|46|). where anomalous response 
first emerges in the frozen state. In this subsection we demonstrate that our dynamical 
equations (jl4|15j) indeed allow for self-consistent stationary state solutions with the 
properties (gTJ). Note that (gZJ) directly imply that ]T r [J + G] _1 ( r ) = i 1 + xY l = 0, 
and hence also a 2 = 0. Insertion of (|4"7j) as ansatze into ()14|15j) now gives: 

= [(1 + G)~ X D{1 + G t )~ 1 G t ](t) - [(1 + G) _1 C](r), (48) 

G(t + 1) = G(t) + a(l + G) _1 (r) + (1 - a)<J(r), (49) 

where S(t) = 1 if r = and 5(r) = otherwise. We define the persistent response 
g = limx^oo T _1 J2 t <t ^( r )' anc ^ sum both sides of from r = to r = f, to get: 

G{£+ 1) = 1 -a + a^(l + G)- 1 (r). (50) 

We conclude that g — 1 — a. From this, in turn, one infers for a < 1 that indeed 
X = J2 T ^( r ) = 00 > which confirms in retrospect that our ansatze (|4*7jl indeed solve our 
dynamical laws (|48|49|) . 



5. The Phase Diagram and its Verification 

In the previous section we have identified three distinct phases, which according to 
the extensive numerical simulations described below exhaust the phase diagram in the 
(a, r)-plane of the present spherical MG. We have also derived explicit expressions for 
the macroscopic order parameters and the volatility in all three phases, and we have 
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Figure 2. Verification of O^AF and O^F transitions. Left pictures: Lagrange 
parameter A, oscillation amplitude c\ and volatility a 2 versus a, for r = 
1.0,0.9,0.8,0.7,0.6,0.5 (from top to bottom in the panels showing c\ and a 2 ). Here 
we expect to sec the O^AF transition. Markers: simulations; solid lines: theory. The 
dashed vertical line marks the predicted location a Ci i of the O— ^AF transition. Right 
pictures: the same observables shown versus a, but now for r = 0.4,0.3,0.2; here we 
expect to see the O— >F transition. Markers: simulations for r = 0.2 (circles), r = 0.3 
(squares) and r = 0.4 (diamonds); solid lines: theory. The lines are continued as 
dashed lines into the opposite phases, where they should no longer be valid. 
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been able to calculate the various phase boundaries in explicit form. Our results may 
be summarised as follows: 

a > max{a c l , a Cj2 (r)} : oscillating phase (O) 

oscillating C(t), finite x 

r < r*, a c>3 (r) < a < a Cj2 (r) : frozen phase (F) 

constant C(t), finite x 

a < min{a Cj i, a Cj 3(r)} : anomalous frozen phase (AF) 



with 



a c ,i = 3 - 2\/2, a Cj2 (r) 



constant C(t), infinite x 



r + l/2rl 2 r 2 



#TT] ' a ' Ar) = r^l- (51) 

The resulting phase diagram is shown in figure Q Let us briefly discuss its main features. 
For r > r* the behaviour of the spherical MG is similar to that of its conventional 
counterparts, with a divergence of x & t some fixed critical a and x remaining infinite 
as a is reduced further to zero. There are crucial differences though: Firstly, in the 
conventional MG persistent oscillations are found only for a < a c , while they decay 
above the transition. The opposite is the case in the spherical model. Secondly, in the 
conventional MG the volatility a 2 is a smooth function of a across the transition. In the 
spherical model we find that the volatility (and the amplitude c\ of the oscillations as 
well) exhibits a jump at a = a c> i for r > r*. The discontinuity of a 2 follows immediately 
from the nonzero value of (j4"Uj) for r > r*, which gives a 2 in the phase O close to the 
O^AF transition, whereas one has a 2 = throughout the AF phase. The magnitude 
of the jump decreases as r is lowered and finally vanishes at r = r*. Below r = r* 
no discontinuities are present and one finds an intermediate regime, where the system 
freezes, but as yet with a finite x- Only as a is lowered further a transition to a frozen 
phase with anomalous integrated response takes place at a = a Cj3 (r), and below a Cj3 (r) 
both the volatility and the normalisation factor A vanish identically. 

We have tested our theoretical predictions against numerical simulations of the 
spherical MG. The data shown in the figures are all obtained from simulations of the 
batch process (|4I5|) with iV = 500 players, and averaged over 20 realisations of the 
disorder (i.e. the realisations of the strategies). All measurements are temporal averages 
over 250 time steps, preceded by 250 'equilibration' steps. We focus on the parameter 
regions where the various phase transitions are predicted to occur and depict the values 
of the stationary order parameters A and c\ as well as the volatility a 2 as indicators for 
the predicted transitions. The precise locations of the various transitions are given in 

Figure El concerns the O— >AF and O— >F transitions. For r > r* w 0.455 we expect 
to see the O— >AF transition. Here our theory predicts that A = a — 1 + 2y / a for 
a > a Cj i, and A = for a < a Ct i. At a Cj ± the volatility and the oscillation amplitude 
C\ should both jump discontinuously to zero. For r < r*, on the other hand we should 
observe the O— >F transition where c\ goes to zero continuously at a = a c ^{r). We 
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Figure 3. Verification of F^AF transition. Lagrange parameter A (left) and volatility 
<7 2 (right) versus a, for r = 0.2 (circles), r = 0.3 (squares) and r = 0.4 (diamonds). 
Solid lines: theory. Both observables are predicted to converge to zero as the transition 
a = a Cl s(r) is approached from above, and to remain zero below the transition. 

also expect A and a 2 to be continuous at this transition, albeit their derivatives with 
respect to a change discontinuously. The data in figure El reveal full agreement between 
theory and simulation (up to finite size effects close to the transitions). Although we 
restrict ourselves to r < 1 in figure El we have verified that the qualitative behaviour of 
the system remains unchanged for larger values of r and that the very good agreement 
between theory and simulation continues to hold for r > 1. Figure H3 concerns the F^AF 
transition, where A and a 2 are predicted to vanish as a approaches d c ^{r) from above. 
Again we find good agreement between theory and numerical experiment. 

6. Conclusions 

In this paper we have introduced a spherical version of the batch minority game, 
with random public information. In this model the non-linear update rule of the 
conventional game is replaced by an iteration prescription which is linear in the 
microscopic degrees of freedom (the point differences <&(£)), complemented by a spherical 
constraint N^ 1 J^qf (£) = 7-2 ■ The spherical MG is designed to be exactly solvable in 
the thermodynamic limit. In terms of the decision making of the individual agents, the 
linearised microscopic dynamical laws corresponds to allowing the agents to play linear 
combinations of their two strategies. The relevant control parameters of the spherical 
MG are the radius r of the sphere to which the dynamics is confined and the ratio 
a = p/N of the number of possible values of the external information over the number 
of agents. 

Using the dynamic mean field theory introduced by De Dominicis we are able to 
perform the average over the disorder and to take the thermodynamic limit. This 
formalism reduces the original iV-agent dynamics to a non-Markovian effective single- 
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agent stochastic process. Like in spherical spin-glass models, the temporal evolution 
of the macroscopic order parameters (the correlation and response functions) can be 
formulated in terms of a pair of coupled iterative equations, without referring to the 
microscopic single effective- agent process. Assuming the existence of a time-translation 
invariant stationary state we are able to solve these equations exactly, and to compute 
the order parameters in the stationary state as well as the stationary volatility at every 
point of the phase diagram without making any approximations. 

We find that, although the update rule is relatively simple compared to the 
conventional MG, the spherical MG displays a remarkably rich structure. Depending 
on r and a the system exhibits three distinct phases, two without anomalous response 
(an oscillating and a frozen state) and a further frozen phase with diverging integrated 
response. As described above the spherical model exhibits some similarities as well as 
intriguing differences compared to conventional MGs. The four main differences are (i) 
the absence of any macroscopic dynamical effect of the choice of the initial microscopic 
state in the spherical game, (ii) the fact that, for any r, the volatility is always zero 
close to a = whereas in the conventional MGs both high-volatility and low-volatility 
solutions can be found, (iii) persistent oscillations in the spherical MG, which increase 
for increasing a and vanish for low a, where in the conventional batch MG persistent 
oscillations can only be found in the low-a regime and (iv) the discontinuous dependence 
of the volatility on a in the spherical MG for r > r*. 

In summary, our study demonstrates that the dynamical rules of the conventional 
MG can be simplified to obtain a completely solvable spherical version, which still 
displays a non-trivial phase diagram and some novel features, which are not observed in 
conventional MGs. It would be interesting to study further the mathematical properties 
of the spherical model, such as the relaxation towards the stationary state and possible 
ageing phenomena. 
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